Localized states influence spin transport in epitaxial graphene 
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We developed a spin transport model for a diffusive channel with coupled localized states that 
result in an effective increase of spin precession frequencies and a reduction of spin relaxation times 
in the system. We apply this model to Hanle spin precession measurements obtained on monolayer 
epitaxial graphene on SiC(OOOl) (MLEG). Combined with newly performed measurements on quasi- 
free-standing monolayer epitaxial graphene on SiC(OOOl) our analysis shows that the different values 
for the diffusion coefficient measured in charge and spin transport measurements in MLEG and the 
high values for the spin relaxation time can be explained by the influence of localized states arising 
from the buffer layer at the interface between the graphene and the SiC surface. 
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The spin dynamics in the diffusive transport regime 
are in general described by the Bloch equation for the 
spin chemical potential fig that describes the three di- 
mensional spin accumulation [T]: 
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with the diffusion coefficient D, the spin relaxation time 
ts and the Larmor frequency cJl = g/iB/fo B, that 
describes the spin precession in a perpendicular mag- 
netic field B with the gyromagnetic factor g (g-factor, 
g = 2 for free electrons) and the Bohr magneton /is- 
Experimentally, spin transport is commonly examined 
by Hanle spin precession measurements (Fig. [TJi) that 
are fitted with the solutions of the Bloch equation 
|l]). Those fits result in D, t$ and the spin relaxation 
length As = Dt$- However, the fits are invariant 
under the transformation D —> cD, t$ —> t~s/c, g — s- eg 
leaving the scaling factor c undefined. To unambigu- 
ously define the parameters, D can be independently 
determined using the diffusion coefficient from charge 
transport measurements Dc and the Einstein relation 
D c = (R sq e 2 i'(E F ))~ 1 [2]. Here R sq is the square 
resistance, e the electron charge and v(Ep) the density 
of states (DOS) of the diffusive channel at the Fermi 
energy. 

Spin transport in graphene has been extensively studied 
in recent years [5lU6|. Due to weak spin-orbit coupling 
g = 2 is commonly assumed to fit Hanle precession 
data (and define c) [3HT5], This was justified for 
exfoliated single layer graphene (eSLG) as it was shown 
that D w Dc 5\. On the contrary, recent results on 
monolayer epitaxial graphene on SiC(OOOl) (MLEG) 
[T71 ITS] show D -C Dc along with very high values for 

T S PUCE]. 

In this letter we introduce a model that can explain the 
apparent difference between D and Dc by the increase of 
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FIG. 1: (Color online) (a) Sketch of the Hanle precession ge- 
ometry with a diffusive channel of width W connected to the 
ferromagnetic spin injector (Fi) and detector (Fd) on distance 
L. The out-of-plane magnetic field B causes the in-plane in- 
jected spins to precess while diffusing through the channel, 
(b) Extension of the Hanle precession geometry with local- 
ized states that are coupled to the diffusive channel. The 
spins can hop into these states and back into the channel 
while the states are not coupled with each other. 



the effective g-factor caused by localized states coupled 
to the spin transport channel. Then we discuss how this 
model reinterprets the results on MLEG from Ref. [TTJ 
and finally, we compare the results on MLEG to Hanle 
precession measurements on quasi-free-standing MLEG 
on SiC(OOOl) (QFMLG) [20]. In this material the 
graphene-like, electrical neutral buffer layer, that is in 
conventional MLEG located between graphene and the 
SiC substrate, is absent [HJ|52]. The presented analysis 
points to the buffer layer as the origin of the localized 
states and by comparing the results from MLEG and 
eSLG we can estimate the spin properties of those states. 
To examine the spin transport properties of graphene, 
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usually the non-local measurement geometry is used, 
consisting of a two dimensional channel with ferromag- 
netic electrodes that inject and detect electron spins in 
the graphene plane [3J [23J (Fig. [iji). To introduce our 
model we extend the description of the spin dynamics 
in this system with localized states in close proximity 
to the channel (Fig. [T|d). The states are electrically 
coupled to the channel and we assume that they are not 
coupled with each other. Therefore spins can hop into 
these states and back into the channel but not directly 
from one localized state into another. 
The spin accumulation in the localized states is rep- 
resented by fi* s and its dynamics can be described by 
a Bloch equation similar to that does not include 
a diffusive term but a term for the coupling to the 
transport channel. 
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In this equation the Lamor precession frequency in the 
localized states uj* l = u* L z = auj^z is introduced, which 
can be different from wl due to a possibly different g- 
factor g* = ag. t s = /3rg is the spin relaxation time 
of the localized states and the term — T((i s — (ig) de- 
scribes the flow of spins from the localized states to the 
transport channel and vice versa with the coupling rate 
r = (RevLs)' 1 , where 1/R is the conductance per unit 
area between the localized states and the graphene chan- 
nel and vls the density of localized states [2"31l2"5] . 
To describe the spin dynamics in the transport channel 
we also have to add on the right side of the Bloch equa- 
tion ([I]) a coupling term in the same way as in equation 
§2§. Therefore we get 
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Here we introduce the factor r/ = VLsl v gr that accounts 
for the different DOS in graphene v gr compared to the 
localized states. 

The two coupled equations ^ and ([3| can be reduced 
to one effective Bloch equation. For this purpose we 

consider the system to be in a steady state with = 0, 
so that equation &2h can be rewritten as (jl* s = a - /Is with 
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As the spin accumulation is purely perpendicular 
to the magnetic field in the Hanle geometry [53] 
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FIG. 2: (Color online) The effective spin relaxation time r| 

(a) and Lamor precession frequency (b) as a function of 
the coupling rate T (black solid lines). The asymptotic values 
in the limit of strong (red, dashed line) and weak coupling 
(blue, dash dotted line). In the shown graphs we keep ts = 
150 ps (gray dotted line in panel (a)), Tg = 5 ns, ujl = 10 GHz 
(for a magnetic field of B as 50 mT, gray dotted line in panel 

(b) ) and r\ — 50 constant. 
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Here we introduce the effective spin relaxation time Tg 
and the effective precession frequency of the system 
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So the effective spin dynamics in the transport channel 
including the localized states can be described by one 
single effective Bloch equation with an effective spin re- 
laxation time and precession frequency. 
The expressions for t s ^ and are plotted in Fig.^as 
a function of the coupling rate V (black solid lines) . Note 
that independent from T the model shows an effective de- 
crease of the spin relaxation time ^ and an effectively 
increased precession frequency ([7]) as the coupled local- 
ized states result in extra relaxation and additional spin 
precession. To further analyze the effective values we can 
consider certain limits. T is inversely proportional to the 
average dwell time in the localized states T oc 1/Tdweii, 
the average time the spins stay in the localized states be- 
fore hopping back into the channel. If we consider weak 
coupling, T <C (Fig. [2j blue dash dotted lines), we 

have long dwell times and therefore Td we ii ^> T s- ^ s 
a consequence all the spins that hop into the localized 
states will relax before returning into the diffusive chan- 
nel and are therefore "lost" for the spin transport. We 
get l/rf f w l/r s + rjT, while oj e L ff = oj l + 0({t* s T) 2 ) 
stays approximately constant. 



If we have strong coupling, r>l /t s (Fig. [2j red dashed 
lines), we have to distinguish between the cases where 
the precession frequency in the localized states is greater 
or lower than the coupling rate. In the case that the 
precession frequency is higher, uj* l ^> T, we get the same 
result for T5 and uj e ^ as for weak coupling. The strong 
precession in the localized states dephases all spins that 
hop into these states so that again all spins hopping into 
these states are lost. The influence on the effective pre- 
cession frequency is therefore marginal. 
The most interesting case is the case of strong coupling, 
T ^> 1/ts, and low precession frequencies, uo* L <C L, 
as this corresponds to the measurements in MLEG (see 
below). We get: l/r e s ff = 1/t s + tj/t s and u e / f = 
ljl + rjij* L . Both values are in this limit independent from 
the coupling rate L (Fig. |2k. 

How does this model relate to the results on spin trans- 
port in MLEG on SiC(OOOl) reported in Ref. [11]? Here an 
increased Tg and a strongly reduced diffusion coefficient 
(D <C Dc) were observed. To understand this obser- 
vation in correspondence to the discussed model, let us 
revisit how the spin transport data is acquired. As men- 
tioned before g = 2 is commonly assumed to fit Hanle 
precession measurements on graphene. In our system 
with an effective precession frequency = £,u>l > 
and hence g e ^ = ^g this assumption presents itself 
wrong. The values that are received by fitting assum- 
ing g = 2 are described by a modified Bloch equation 
that we receive by dividing ^ by the scaling factor £: 

= £ mod VVs "4 + ^ x Ms (8) 
T s 

with D mod = D/£ and T™ od = ^T e J s ' . The effective spin 
relaxation time of the system including the localized 
states can be obtained by either assuming D = Dc 
for the fit or assuming g — 2 and correcting the spin 
relaxation time with t% = T s nod /£,. The enhanced 
value measured for ts in Ref. [TT] is therefore not an 
intrinsic property of MLEG on SiC(OOOl) but is based 
on assuming a value for g of graphene without taking 
the influence of the localized states into account. 
Note that the measured spin relaxation length does 
not change when assuming a different g-factor as 

X mod = ( D mod T mody/2 = _ A |?/_ jjg^ 

the spin relaxation length in the system is reduced by 
the influence of the localized states. 
To show the effect of a modified g-factor let us revisit 
the data in Ref. [TT] The narrow Hanle precession 
measurement in the center of Fig. [3] is performed on 
MLEG on SiC(OOOl) at room temperature (RT). The 
fit to this data (assuming g — 2) gives ts = 1.3 ns 
and D = 2.4 cm 2 /s, resulting in As = 0.56 fj,m. This 
fit shows the increased value of ts compared to values 
obtained on eSLG [5] [TT] but also the strongly reduced 
values for D compared to the value of Dc ~ 190 cm 2 /s 
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FIG. 3: (Color online) Hanle precession measurements per- 
formed at RT with parallel (tt) and antiparallel alignment 
(t4) °f the inner electrodes. The measurements with the nar- 
row curve in the middle (light colors, scale on the left axis) 
is taken from Ref. [TT] and were performed on MLEG on a 
L — 1.2 jj,m long and W = 0.7 fim wide strip. The broader 
measurements enclosing the MLEG measurements (dark col- 
ors, scale on the right axis) were performed on QFMLG with 
L = 1.5 fim and W = 1 /im. The fits to the solutions of the 
Bloch equation |T| are plotted in gray. For both sets of mea- 
surements a constant background resistance was subtracted. 

obtained in charge transport measurements on the same 
sample and compared to D ~ 200 cm 2 /s typically 
measured on eSLG [5j. 

To find out where the predicted localized states originate 
from, we prepared and measured spin transport samples 
on quasi- free-standing MLEG (QFMLG). To obtain 
QFMLG only the electrical neutral buffer layer and 
no graphene layers are grown on SiC(0001) and the 
sample is then intercalated by hydrogen as described in 
Refs. dTJ and Thus the buffer layer decouples from 
the substrate and converts into a graphene layer while 
dangling bonds of the SiC surface are passivated. We are 
therefore left with a single graphene layer directly on the 
SiC(0001) surface. The experiments were performed in 
the same way as described for MLEG samples in Ref. [TT] 
Fig. [3] shows a Hanle precession curve measured on a 
QFMLG strip at RT next to the measurement on MLEG 
from Ref. 11 This Hanle curve shows clearly a big 
difference compared to the measurement on MLEG on 
SiC(0001) in the same figure. The non-local resistance 
[3] [TT] changes slower with the magnetic field [2"o] . 
comparable to measurements performed on eSLG [5J. 
The fit (assuming g = 2) gives rg = 33.6 ± 0.9 ps and 
D = 75 ± 2 cm 2 /s and therefore A s = 0.50 ± 0.01 //m. 
These values are similar to a low quality eSLG sample 
as ts is reduced by about a factor 4 and D by a factor 
of approximately 3 compared to typical values obtained 
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on eSLG [5]. Compared to MLEG we see an increase 
of D by a factor of ~ 30 and a decrease of ts by 
about 40 times. And in contrast to the data obtained 
on MLEG we see in charge transport measurements 
on similar QFMLG samples a diffusion coefficient of 
D c « 45 cm 2 /s ~ D [37]. To obtain D c we use the 
square resistance of R sq pa 3.5 kfi and a hole charge 
carrier density of p pa 6 x 10 12 cm~ 2 obtained in Hall 
measurements consistent with results from Ref. [22] 
Comparing the results on the two types of graphene on 
SiC(OOOl) it is interesting to see the striking difference 
of the spin relaxation times and diffusion coefficients 
obtained in spin transport measurements but even more 
important that we see Dc pa D in QFMLG as we would 
expect for graphene [5J. This points to the fact that 
here is no effect of the localized states. Therefore, the 
localized states have their origin in the interface between 
the graphene layer and the SiC substrate as this is 
the only structural property that is altered between 
conventional MLEG and QFMLG. Hence, the states 
could be in the dangling bonds or in the buffer layer. 
The strong difference of D vs Dc reported in Ref. [TTJ 
points to a strong change in ujl- Within our model, this 
is the case for strong coupling of the localized states to 
the channel (see eq. ^ and Fig. ^p) . The coupling is 
strongest if the localized states are located in the buffer 
layer as this one is closest to the channel. If we assume 
that the coupling of the states in the buffer layer to the 
channel is comparable to the coupling between adjacent 
layers in graphite and considering rj ~ 50 (see below) 
we get r ~ 2 x 10 13 s _1 [25]. For this value the strong 
coupling limit seems justified (see Fig.[2|. 
Now we can evaluate the model and characterize the 
localized states by comparing the fitting results on 
MLEG on SiC(OOOl) from Ref. [TTJ with data obtained 
on other types of monolayer graphene. To compensate 
for different D^-values obtained in charge transport 
measurements on QFMLG and conventional MLEG, we 
use data on eSLG [3] to compare with the fitting results 
on conventional MLEG [TT]. In the limit of strong 
coupling we get based on £ = 1 + arj. We assume 
that the g-factor in the localized states is equal to the 
graphene channel (a — 1). Using the typical eSLG 
values, ts = 150 ps and Dc = D = 200 cm 2 /s, as the 
graphene values in the absence of localized states and 
the values obtained on MLEG as D mod and r s nod we get 
£ pa n pa 80 at RT. With this result and r^ orf /r s pa 9 we 
obtain t s /ts — (3 ~ 10. Hence, at RT spins relax in the 
localized states with t s pa 1.5 ns about 10 times slower 
than in the graphene channel. This enhanced value is 
very reasonable for a confined state in a material with 
low spin-orbit coupling. 

The presence of localized states can also explain the 
temperature dependence of the spin transport properties 
in Ref. [TT] in contrast to only negligible change for 
eSLG [3]. By assuming the same values as before as the 



values of D and t$ in the absence of localized states, 
we get with the data for MLEG at 4 K r\ pa 45 and 
(3 sa 22 (rg sa 3.3 ns). These results imply that at low 
temperature there are less localized states accessible 
and those states have a longer spin relaxation time. 
By assuming a Boltzmann distribution we get from 
the change in t) an activation energy for the states of 
E a pa 15 meV. 

Within our analysis, 77 describes the ratio of the DOS 
in the localized states and the graphene channel. With 
r) up to 80 we need a high density of localized states in 
our system. In MLEG we have with a electron charge 
carrier density of n « 3 x 10 12 cm -2 [TT] a DOS of 
10 13 eV 'cm 2 . With a density of carbon 
atoms in the graphene-like buffer layer of 3.8 x 10 15 cm~ 2 
and assuming that every carbon atom contributes one 
localized state, we get r\ — 80 if these states are e.g. 
uniformly distributed over an energy range of pa 1 eV. 
Those localized states can be the origin of the strong 
doping observed in MLEG on SiC(0001) [25]. 
The observed increase of g in MLEG could in prin- 
ciple also be related to magnetic moments induced 
by the buffer layer or dangling bonds on the surface 
of SiC (0001) as described for hydrogenated graphene 
in Ref. [14] We argue that this does not apply here 
since: i) The effect in MLEG is stronger at RT than at 
4 K while the effect in hydrogenated graphene is only 
observed at low temperature, ii) We do not see any 
effects resulting from randomized magnetic moments 
at low magnetic fields like the "dip" in the spin-valve 
measurements in Ref. O iii) The increase of g in MLEG 
is much bigger than in hydrogenated graphene. 
To summarize, we developed a spin transport model 
for a diffusive channel with coupled localized states 
that results in an increased effective g-factor and a 
reduced spin relaxation time for the transported spins. 
This model reinterprets the data from Ref. [TTJ where 
an enhanced spin relaxation time and a reduced spin 
diffusion coefficient were observed. By comparing the 
data from Ref. [TTJ to new measurements on QFMLG 
and typical values on eSLG we could identify the 
buffer layer as possible source for the localized states 
and the measurements can be related to a g-factor 
of g e ff = (45 — 80)g. Finally we use the model to 
characterize the spin properties of the localized states in 
the buffer layer of MLEG on SiC(0001). 
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SUPPLEMENTAL MATERIAL FOR 
"LOCALIZED STATES INFLUENCE SPIN 
TRANSPORT IN EPITAXIAL GRAPHENE" 

Coupling rate between the localized states and 
the graphene channel 

To estimate the coupling rate T between the local- 
ized states and the graphene channel, we can set up a 
simplified model based on the coupling between adja- 
cent graphene layers in graphite, as these layers have the 
same or similar physical distance as the buffer layer to 
the graphene layer. 

In graphite, the conductance in z-direction perpendicular 
to the layers is per layer an, = a gr /{C,d) where a gr is the 
in-plane conductivity of a graphene layer, d the distance 
between two layers (or between the graphene layer and 
the localized states) and £ « 100 the ratio between the 
conductivity within the layers and perpendicular to them 



Using the Einstein relation with v gr the DOS and D the 
diffusion coefficient of the graphene channel we get: 



y Vg^P_l = dV 
vls d 2 C dt 



(10) 



This equation includes the ratio of the DOS of the local- 
ized states and the graphene channel r\ = v^s/vgr that 
was discussed in the main text of this letter. 



With the coupling rate r ~ X ^nr we receive 
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r ~ v d? c 



(11) 



With this model we get for bilayer graphene with £ = 
100, v gr = vlSi d = 0.3 nm and the typical graphene 



value D « 0.02 m 2 /s 
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10 15 s" 1 . 
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We can now calculate for a current In in z-direction: 
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dQ 
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Here V — fi/e is the voltage between the localized states 
and the channel, proportional to the difference in the 
chemical potential, A the area through which the current 
flows, Q is the total charge that flows, N the number of 
charge carriers, d is the distance to and v^s the density 
of states (DOS) of the localized states, and e the electron 
charge. 



For our system we have r/ ~ 50 while the other parame- 
ters stay the same and get therefore 



T LS «2x 10 13 s" 1 . 



(13) 



This value gives the order of magnitude of the coupling 
rate between the localized states and the graphene 
channel. With this value we are in the limit of strong 
coupling of the system as depicted in Fig. 2 of the main 
text. 



